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Abstract 

Based on a new approximation method, namely pseudospectral method, a solution for the 
three order nonlinear ordinary differential laminar boundary layer Falkner-Skan equation 
has been obtained on the semi-infinite domain. The proposed approach is equipped by the 
orthogonal Hermite functions that have perfect properties to achieve this goal. This method 
solves the problem on the semi-infinite domain without truncating it to a finite domain and 
transforming domain of the problem to a finite domain. In addition, this method reduces 
solution of the problem to solution of a system of algebraic equations. We also present 
the comparison of this work with numerical results and show that the present method is 
applicable. 
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1. Introduction 

There are many problems in science and engineering arising in unbounded domains. 
Spectral methods are famous ways to solve these kinds of problems. The most common 
approach on spectral methods, that is used in this paper too, is through the use of polyno- 
mials/functions that are orthogonal over unbounded domains, such as the Hermite and the 
Laguerre polynomials / functions [lMs|] . 

Guo |9l-lll| proposed a method that proceeds by mapping the original problem in an 
unbounded domain to a problem in a bounded domain, and then using Jacobi polynomials 
to approximate the resulting problems. 

Another method is replacing infinite domain with [-L, L] and semi-infinite interval with 



[0, L] by choosing L, sufficiently large. This method is called domain truncation 12 
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Using rational approximation is another effective direct approach for solving such prob- 
lems. Christov 13[ and Boyd 14[ ll5| developed some spectral methods on unbounded 
intervals by using mutually orthogonal systems of rational functions. Boyd [l5[ defined a 
new spectral basis, named rational Chebyshev functions on the semi-infinite interval, by 
mapping to the Chebyshev polynomials. Guo et al. [l6[ introduced a new set of rational 
Legendre functions which are mutually orthogonal in L 2 (0, +00). They applied a spectral 
scheme using the rational Legendre functions for solving the Korteweg-de Vries equation on 



the half line. Boyd et al. [17] applied pseudospectral methods on a semi-infinite interval 



and compared rational Chebyshev, Laguerre and mapped Fourier sine. 
Parand et al. 



18H23J applied spectral method to solve nonlinear ordinary differential 



equations on semi-infinite intervals. Their approach was based on rational Tau and pseu- 
dospectral methods. 



The concept of boundary layer flows [2J| of an incompressible fluid has several engineer- 
ing applications, such as aerodynamic extrusion of plastic sheets and cooling of a metallic 
plate in a cooling bath, cooling and heating of the plate by free convection current [2o| . 
Such flows are also encountered in the glass and polymer industries. The effect of ther- 
mal radiation absorption on an unsteady free convective flow past a verti cal p late has been 
studied in the presence of a magnetic field and constant wall heat flux by 25( also Makinde 
j26[ investigated numerically the hydromagnetic mixed convection flow of an incompress- 
ible viscous electrically conducting fluid and mass transfer over a vertical porous plate with 
constant heat flux embedded in a porous medium by Newton-Raphson shooting method 
along with fourth-order Runge-Kutta integration algorithm. The case of mixed-convective 
boundary layer flow past a vertical porous plate embedded in a saturated porous medium 
with a constant heat flux and mass transfer in the presence of a magnetic field and with 
heat absorption has been studied by 27 . 

The steady laminar flow passing a fixed wedge was first analyzed in the early 1930s by 
Falkner and Skan 28) to illustrate the application of Prandtl's boundary layer theory. The 
Magnetohydrodynamic (MHD) flows as boundary layer flows have a great significance both 
from a mathematical as well as a physical standpoint. Such flows are very important in 
electromagnetic propulsion. The MHD systems are used effectively in many applications 
including power generators, pumps, accelerators, electrostatic filters, droplet filters, the 
design of heat exchangers, the cooling of reactors etc. 
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1.1. MHD flows 

Consider the boundary layer flow of an electrically conducting viscous fluid. A magnetic 
field B(x) acts transversely to the flow. The induced magnetic field is neglected by choosing 
small magnetic Reynolds number assumption. Furthermore, the electric field is absent. The 
relevant problem is 



du dv q 
dx dy 



(1) 



du du y.du d 2 v 
dx dy dx dy 2 



0B 2 



[u 



U), 



(2) 



u = 0, v = at y = 0, w = t/(x), as y —>■ oo, (3) 

where Q 



E/"(z) = ax m (4) 
and 

= 5 x (m - 1)/2 , (5) 

in which it and u are the velocity components, U is the inherent characteristic velocity, v is 
a kinematic viscosity, a is the electrical conductivity, p is the fluid density, B and £?o are 
the magnetic field and externally imposed magnetic field in the y-direction respectively. 
Defining 



lm + 1 I U 

T 




—y, f = ,— vW/(r), (6) 
vx V m + 1 



V 2 V x m + 1 

the continuity equation is identically satisfied and Eq. |5] and boundary conditions [3] reduce 
to the following form 

+ /0 + fl[l-(|) 2 ]-M 2 (^-D=O, (8) 

/(0) = /'(0) = 0, /'(+oo) = l, (9) 

where (3 = and M 2 = 2aB^/ pa(l + m). Now our interest is to find solution of Eq. (jHJ) 
for wedge in the accelerated flow (m > 0, (3 > 0) and decelerated flow (m < 0, /3 < 0) with 
separation. 



MHD boundary layer flows have been studied by several researchers. Yih [3JJ and Ishak 
et al. [32| transformed the partial differential boundary layer equations into the nonsimilar 
boundary layer equations and a system of ordinary differential equations respectively, then 



they used Keller box method to solve them. Hayat et al. [33j , Rashidi [34| and Abbasbandy 



et al. [35|,[36( solved MHD boundary layer flow by modified Adomian decomposition method, 
DTM-Pade, Homotopy analysis method and Hankel-Pade respectively. 

The purpose of this paper is to develop a pseudospectral method that is more direct and 
simpler than the other spectral methods. The method of this paper differs from previous 
methods that have been used because it approximates the solution directly. Thus, it does 
not truncate the semi-infinite domain to a finite domain, does not introduce the unknown 
finite boundary, does not impose an asymptotic condition at this unknown boundary, and 
does not use transformation to map the physical domain. 
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This paper is organized as follows: 
In section El the properties of Hermite functions and the way to construct the pseudospectral 
approach for this type equation are described. In section [3] the proposed method is applied 
to solve MHD Falkner-Skan equation, and a comparison is made with existing numerical 
solutions that were reported in other researches, and we have a discussion in the final 
section. 



2. Hermite functions pseudospectral method 

Spectral methods have been successfully applied in the approximation of differential 
boundary value problems defined in unbounded domains. For problems whose solutions 
are sufficiently smooth, they exhibit exponential rates of convergence/spectral accuracy. 
We can apply different spectral methods that are used to solve problems in semi-infinite 
domains. One of these approaches is using Laguerre and Hermite polynomials/functions 
@, S^S, 37 39] . Guo jl] suggested a Laguerre- Galerkin method for the Burgers equation and 



Benjamin-Bona-Mahony (BBM) equation on a semi- infinite interval. In [7j Shen proposed 
spectral methods using Laguerre functions and analyzed for model elliptic equations on reg- 
ular unbounded domains. Siyyam Q applied two numerical methods for solving initial value 
problem differential equations using the Laguerre Tau method. Maday et al. |6[ proposed a 
Laguerre type spectral method for solving partial differential equations. Funaro and Kavian 
[2| considered some algorithms by using Hermite functions. Recently Guo [H developed 
the spectral method by using Hermite polynomials. However, it is not easy to perform the 
quadratures in unbounded domains, which are used in the Hermite spectral approximations. 
So the Hermite pseudospectral method is more preferable in actual calculations. Guo [39I ] 
developed the Hermite pseudospectral method for the Burgers equation on the whole line. 



Guo et al. [37J] considered spectral and pseudospectral approximations using Hermite func 



tions for PDEs on the whole line to approximate the Dirac equation. Bao and Shen [38 
proposed a generalized-Laguerre-Hermite pseudospectral method for computing symmetric 
and central vortex states in Bose-Einstein condensates (BECs) in three dimensions with 
cylindrical symmetry. 

Pseudospectral methods has become increasingly popular for solving differential equa- 
tions and also it is very useful in providing highly accurate solutions to differential equations. 
In this paper, we employ the Hermite functions pseudospectral method which we denote HFP 
to solve MHD Falkner-Skan initial value problem directly. 

2.1. Properties of Hermite functions 

In this section, we detail the properties of the Hermite functions that will be used to 
construct the HFP method. First we note that the Hermite polynomials are generally not 



suitable in practice due to their wild asymptotic behavior at infinities |40|, |41 . 

Hermite polynomials H n (x), n > 0, are the eigenfunctions of the singular Sturm-Liouville 

problem in 

Hl(x) - 2xH' n (x) + 2nH n (x) = (10) 
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Hermite polynomials with large n can be written in direct formula as follow: 



tt , x T(n+1) r 2 /2 . r- nir. 



~ n n/2 e x2/2 cos( v / 2^TTa: 



nir , 



Hence, we shall consider the so called Hermite functions. The normalized Hermite functions 
of degree n are defined by 



HJx) 



-.e' x2/2 H n (x), n>0,xe 



V2 n n! 

Clearly, {H n } is an orthogonal system in L 2 (M), i.e., 



r 

J — oo 



H n (x)H m (x)dx = ^8 ri 



where 5 nm is the Kronecker delta function. 

In contrast to the Hermite polynomials, the Hermite functions are well behaved with the 
decay property: 



\H n (x) 



0, as \x\ 



oo, 



and the asymptotic formula with large n is 



H n (x) ~ n~4 cos(v / 2rT+Ta; — ) 

The three-term recurrence relation of Hermite polynomials implies 



H n+1 (x) = x^j ^H n (x) - y/^Hn-iix), n > 1, 
H (x) = e~ x2 l 2 , H^x) = ^/2xe~ x2 l 2 . 
Using recurrence relation of Hermite polynomials and the above formula leads to 

H'Jx) = V2^H n ^(x) - xH n (x) 



and this implies 



H' n (x)H' m (x)dx 



— - — j > m — n — 2, 

0r(n + |), m = n, 

_ V(n + l)(n + 2)^ m = n + 2? 

0, otherwise. 
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Let us define 

P N ■= { u : u = e~ x2/2 v,\/v G P N }. 

where P/v is the set of all Hermite polynomials of degree at most N. 
We now introduce the Gauss quadrature associated with the Hermite functions approach. 
\ N 

■ n 



Let {xj}jL be the Hermite-Gauss nodes and define the weights 



Wj = , < j < N. 



Then we have 

N 



3=0 



For a more detailed discussion of these early developments, see the |42|, |43 



2. 2. Approximations by Hermite functions 
Let us define A := {x\ — oo < x < oo} and 

H N = span{H (x), H x (x), H N (x)} 

The L 2 (A)-orthogonal projection ^ : L 2 (A) — > is a mapping in a way that for any 
v G /. 2 (.\). 

< £nV — v , 4> >= 0, V0 G Hn 
or equivalent ly, 

N 

£nv(x) = 22viHi(x). 

1=0 

To obtain the convergence rate of Hermite functions we define the space H r A (N) defined by 

H A (A) = {v\v is measurable on A and \\v \\ rA < oo}, 

and equipped with the norm ||i>|| r)J 4 = \\A r v \\. For any r > 0, the space H A (A) and its norm 
are defined by space interpolation. By induction, for any non-negative integer r, 

r 

A r v(x) = ^(x 2 + l) {r - k)/2 Pk (x)d k x v(x), 

k=0 

where Pk{x) are certain rational functions which are bounded uniformly on A. Thus, 

. 1/2 

\v\\ r ,A < C I J2 II + l) ir - k)/2 Pk(x)d k x V 
\k=0 

Theorem : For any v G H A (A), r > 1 and < /i < r, 

|| Znv-vW^ cN^-W II v || r)A . (11) 

Proof. A complete proof is given by Guo et al. [37|. also same theorem has been proved by 
Shen et al. 40 . 



2. 3. Hermite functions transform 

As mentioned before, Falkner-Skan equation is defined on the interval (0, +oo); but we 
know properties of Hermite functions are derived in the infinite domain (— oo, +oo). 
Also we know approximations can be constructed for infinite, semi-infinite and finite inter- 
vals. One of the approaches to construct approximations on the interval (0, +oo) which is 
used in this paper, is to use a mapping, that is a change of variable of the form 

u = <P(z) = ^\n(z). (12) 

where k is a constant. 

The basis functions on (0, +oo) are taken to be the transformed Hermite functions, 

H n (x)=H n (x)o<j ) (x) = H n (<f ) (x)). 
where H n (x) o 0(x) is defined by H n ((p(x)). The inverse map of u = <fi(z) is 

z = = e kul . (13) 

Thus we may define the inverse images of the spaced nodes as 

T = {0 _1 (t) : -oo < t < +00} = (0, +00) 

and 

x j = (f>-\x j ) = j = 0,1,2,... 

Let w(x) denote a non-negative, integrable, real-valued function over the interval T. We 
define 

l4(r) = {v : T — > R I v is measurable and || v\\ w < 00} 

where 

1 

00 \ 5 







v(x) I w(x)dx 

is the norm induced by the inner product of the space L^(r), 

POO 

<u,v> w = / u{x)v (x)w(x)dx. (14) 







Thus {H n (x)} n( zfq denotes a system which is mutually orthogonal under ffT4l) . i.e., 

< H n (x),H m (x) > W ( X )= V7r5 nmj 

where w(x) = 1/x and 5 nm is the Kronecker delta function. This system is complete in 
L^(r). For any function / G L^(T) the following expansion holds 



N 
k=0 

with 

< f(x),H k (x) > w{x) 



fi 



k 



H k (x) 



I 2 

\w(x) 



Now we can define an orthogonal projection based on transformed Hermite functions as 

below: 

Let 

K N = span{H (x), H^x), H n (x)} 
The L 2 (r)-orthogonal projection : L 2 (T) — > "Hat is a mapping in a way that for any 

2/a 2 (r), 

< £nV -y,(/>>= 0, V0 e H N 
or equivalently, 

TV 

£isry{x) =y]aiHi(x). (15) 

i=0 

2. 4- Domain scaling 



It has already been mentioned in [44j that when using a spectral approach on the whole 
real line M. one can possibly increase the accuracy of the computation by a suitable scaling 
of the underlying time variable t. For example, if y denotes a solution of the ordinary 
differential equation, then the rescaled function is y(t) = y(j), where I is constant. 
Domain scaling is used in several of the applications presented in the next section. For more 
detail we refer the reader to 45 . 



3. Solving Falkner-Skan wedge flow equation with Hermite functions 

Here, we solve Falkner-Skan wedge flow equation for m = —3/5 with M = 5, 10, 15,20 
and m = 2 with M = 5, 10, 50, 100 by pseudospectral method. 

We note that the Hermite functions are not differentiable at the point x = 0, therefore 
to approximate the solution of Eq. (jSJ) with the boundary conditions Eq. §§§ we multiply 
operator Eq. f|T5l) by x 2 /(x + 1) to satisfy /(0) = /'(0) = in Eq. (jHJ) and also we construct 
a function p(x) that satisfy /'(+oo) = 1 in Eq. ([9]). This function is given by 



s 



(16) 




where A is constant to be determined. 

Therefore, the approximate solution of f(x), in Eq. (jHJ) with boundary conditions Eq. (jSJ 
is represented by 

InJ(x) = p(x) + — —£ N f(x), 

X +1 

To apply the pseudospectral scheme, we construct the residual function by substituting f(x) 
by £zv/(zO m the Falkner-Skan Eq. (jSJ): 

Resi(x) = + (£ N f{x/l) ) I I (17) 

d&/W|) )2 \ _ M2 / d6v/(x/Z) _ 
dx I I dx 

where / is a constant that is defined in domain scaling description, before. £nJ{x) will be a 
good approximation of function f(x), if the residual function is zero on the whole domain. 
In other words, we should select coefficients a«s so that the residual function tends to zero 
on most of the domain. The pseudospectral scheme for Falkner-Skan equation is to find 
CnI{x) such that 

Res l (x j ) = 0, j = 0,...,N + l, (19) 

where the XjS are N + 2 transformed Hermite-Gauss nodes by Eq. [131 This generates a set 
of N + 2 nonlinear equations that can be solved by Newton method for unknown coefficients 
OjS and A. 

The physical quantities of interest which is represented by value of f"(0) is the skin 
friction coefficient. The approximations of the f"(0) obtained by this method and numerical 
value for m = -3/5 with M = 5, 10, 15, 20 and m = 2 with M = 5, 10, 50, 100 are listed 
in Tables [TJ and [2] respectively. 

The velocity boundary layer of the wedge with m = —3/5 and m = 2 for various various 
M are shown in figures [1] and [2] respectively. The value of f"(0) can be seen to increase with 
M. Hence, the presence of a magnetic field also increases the skin friction. 

Logarithmic graphs of absolute coefficients |a»| of Hermite function in the approximate 
solutions for m = —3/5 with M = 50 and m = 2 with same M are shown in Figures [3] and 
H] respectively. The graphs illustrate that the method has an appropriate convergence rate. 



4. Conclusion 

The FalknerSkan non-linear ordinary differential equation arises in the study of laminar 
boundary layers. This problem is the subject of an extensive research to solve because of 
its importance in boundary layer theory, a solution for the non-linear ordinary differential 
laminar boundary layer Falkner-Skan equation has been obtained by pseudospectral method 
on semi-infinite domain, also in this work we concentrated on the boundary conditions of 
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/(O) = /'(0) = 0,/'(+oo) = 1, which is corresponding to a fixed and impermeable wedge 
flow. The physical quantities of interest which is represented by the value of f"(0) is the 
skin friction coefficient. According to the tables and figures value of /"(0) can be seen to 
increase with M. Hence, the presence of a magnetic field also increases the skin friction. 

Most numerical methods reported in the literature thus far are based on transforming 
maps of the physical domain [0, +oo) to the finite domains, shooting methods or finite- 
difference methods obtained by first truncating the semi-infinite physical domain of the 
problem to a finite domain at an unknown finite boundary, which is determined as part of 
the solution by imposing an "asymptotic boundary condition" at this boundary. 

The method presented in this paper used a set of Hermite functions and solved this 
problem on the semi-infinite domain without truncating it to a finite domain, imposing the 
asymptotic condition transforming and transforming domain of the problem. These func- 
tions are proposed to provide an effective but simple way to improve the convergence of the 
solution by pseudospectral method. The validity of the method is based on the assump- 
tion that it converges by increasing the number of Gauss points. Through the comparisons 
among the numerical solutions of Asaithambi and the current work, it has been shown that 
the present work has provided acceptable approach for Falkner-Skan equation; also it was 
confirmed by the theorem and logarithmic figures of absolute coefficients that this approach 
has exponentially convergence rate. In total, an important concern of spectral methods is 
the choice of basis functions; the basis functions have three number of properties: easy to 
computation, rapid convergence and completeness, which means that any solution can be 
represented to arbitrarily high accuracy by taking the truncation N sufficiently large. 
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Tabic 1: 


The comparison of HFP approximate solution of /"(0) with numerical results ]M 


| when m=-3/5 


M 


N 


k 


I 


Present method 


Numerical 


5 


20 


2 


1.658 


4.60075494 


4.60075494 


10 


15 


2 


1.296 


9.80646420 


9.80646420 


15 


15 


1 


1.089 


14.87167484 


14.87167484 


20 


15 


1 


1 


19.90393701 


19.90393701 


50 


20 


2 


1.336 


49.96165233 


49.96165233 



Table 2: 


The comparison of HFP approximate solution 


of /"(0) with numerical results [46| when m = 2 


M 


N 


k 


/ 


Present method 


Numerical 


5 


30 


3 


1.194 


5.19095945 


5.19095945 


10 


30 


2 


1.112 


10.09677545 


10.09677545 


50 


30 


2 


0.904 


50.01944071 


50.01944071 


100 


30 


2 


0.616 


100.00972170 


100.00972170 
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Figure 1: The velocity boundary layer of the wedge with m = —3/5 and various M 
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Figure 2: The velocity boundary layer of the wedge with m — 2 and various M 
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Figure 3: Logarithmic graph of absolute coefficients |aj| of Hermite functions in the approximate solution 
for to = -3/5 and M = 50 




Figure 4: Logarithmic graph of absolute coefficients \ai\ of Hermite functions in the approximate solution 
for to = 2 and M = 50 
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